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We identify a new universality class in one-dimensional driven open quantum systems with a dark 
state. Salient features are the persistence of both the microscopic non-equilibrium conditions as 
well as the quantum coherence of dynamics close to criticality. This provides a non-equilibrium 
analogue of quantum criticality, and is sharply distinct from more generic driven systems, where 
both effective thermalization as well as asymptotic decoherence ensue, paralleling classical dynamical 
criticality. We quantify universality by computing the full set of independent critical exponents 
within a functional renormalization group approach. 


Introduction There has been a surge of activity in a 
broad spectrum of experimental platforms, which imple¬ 
ment driven open quantum systems. In such systems, co¬ 
herent and driven-dissipative dynamics occur on an equal 
footing. While such a situation is reminiscent of conven¬ 
tional quantum optics, the systems in point are set apart 
from more traditional realizations by a large, continuous 
number of spatial degrees of freedom, giving rise to gen¬ 
uine driven many-body systems. Indeed, experimental 
realizations range from exciton-polariton systems [1, 2] 
over ultracold atoms [3-5] , large systems of trapped ions 
[6, 7] and photon Bose-Einstein-Condensates [8] to mi¬ 
crocavity arrays [9, 10]. The driven nature at the mi¬ 
croscale leads to an intrinsic non-equilibrium (NEQ) sit¬ 
uation even in stationary state due to the explicit break¬ 
ing of detailed balance. At this microscopic level, the 
dynamics of such systems is Markovian, i.e. memoryless 
in time. This is not in fundamental contradiction to gen¬ 
uine quantum effects playing a role, as has been demon¬ 
strated theoretically [11, 12] and experimentally [13, 14] 
in many-body systems, where phase coherence or entan¬ 
glement ensue in the stationary state of tailored driven- 
dissipative evolution. For the universal critical behavior 
of such systems, however, despite the fact that they are 
“made of quantum ingredients”, the Markovian character 
generically leads to a NEQ analogue of classical dynami¬ 
cal criticality: typically, as the result of dissipation, phase 
transitions in driven-dissipative systems are governed by 
an emergent effective temperature together with the loss 
of quantum coherence, and their bulk critical behaviour 
is captured by equilibrium universality classes [15-24]. 

This sparks a natural and fundamental set of ques¬ 
tions: Given the intrinsic quantum origin, together with 
the flexibility in designing such systems, to which ex¬ 
tent can effects of quantum mechanical coherence persist 
asymptotically at the largest distances in the vicinity of 
a critical point? And if so, what are the precise par¬ 
allels and differences to criticality in closed equilibrium 
systems at zero temperature? In other words, is there a 
driven analogue of quantum critical behavior? 

In this work, we address these questions driving a 
one dimensional open Bose gas with a strong Markovian 


quantum diffusion, implemented, e.g., with microcavity 
arrays. When diffusion dominates over the Markovian 
noise level induced by the environment, a novel critical 
regime associated to NEQ condensation can be realised, 
where the coherent quantum mechanical origin of the 
system and the NEQ driven nature persist at infrared 
scales. Our results establish a new driven quantum uni¬ 
versality class in one dimension, which we characterize 
by computing the full set of static and dynamical crit¬ 
ical exponents. In particular, we obtain the following 
key results: (i) New non-equilibrium fixed point. The 
fixed point (FP) associated to quantum NEQ conden¬ 
sation cannot be mapped to the classical FP of driven- 
dissipative condensation - as evidenced from novel scal¬ 
ing of the correlation length close to criticality, and we 
identify a scaling regime where it governs the fluctuation 
dominated renormalization group (RG) flow. Due to the 
fine-tuning of the Markovian noise level in addition to 
the mass gap, it is less stable - in an RG sense - than a 
classical fixed point, in analogy to the double fine-tuning 
of mass and temperature to zero necessary to reach an 
equilibrium quantum critical point [25]. (ii) Absence of 
decoherence. In classical equilibrium and driven critical 
behavior, decoherence causes the fadeout of all coherent 
couplings in the infrared RG flow, leading to a FP where 
only purely dissipative dynamics persists. In contrast, 
the anomalous dispersion relation of the critical modes 
proper of the novel NEQ fixed point, manifests the simul¬ 
taneous presence of coherent and diffusive processes. The 
absence of decoherence is reflected in the degeneracy of 
the two critical exponents encoding kinetic mechanisms, 
(iii) Absence of asymptotic thermalization. Many driven 
systems exhibit effective thermal behavior at low frequen¬ 
cies [15-22, 24]. The present system does not show this 
property, and this is characterized by the non-thermal 
character of the distribution function, as well as univer¬ 
sally by a new independent critical exponent entering the 
NEQ fluctuation-dissipation relation. A hallmark of the 
interplay of these effects is an oscillatory behaviour of the 
spectral density as a consequence of the (iv) RG limit- 
cycle behavior of the complex wave function renormal¬ 
ization coefficient. 
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FIG. 1. (Color online) A one-dimensional array of microwave 
resonators, coupled to an array of superconducting qubits 
(blue dots), which can decay with rate y (cavity bosons can 
tunnel among neighbouring sites - yellow arrows). Each pair 
of adjacent photonic modes interact with a single qubit via the 
dipole term, TL ~ flcrL (bi — bi+ 1) + h.c.\ bi are the bosonic an¬ 
nihilation operators for the cavity modes and the local qubit 
Hamiltonians are proportional to the cr| Pauli matrix. For 
an energy scale separation r/ Q, the qubit dynamics can 
be adiabatically eliminated [37]. This gives rise to Lindblad 
operators proportional to ~ bi — K+i, which in the continuum 
limit yields Ld(x). It imprints an additional diffusion on the 
propagation of bosons and, crucially for this work, gives rise 
to a scaling of noise level ~ q 2 , as discussed after Eq. (2). 


Finally, notice that the noise in our system is Marko¬ 
vian, in contrast to previous realizations of NEQ quan¬ 
tum criticality [17, 18], while the quantum nature of the 
novel critical regime sharply sets apart our scenario from 
other NEQ fixed points (FP), as occurring in surface 
growth [26] , directed percolation [27] , or turbulence [28- 

31]- 

The results presented here are obtained within a func¬ 
tional renormalization group approach [32-34] based on 
the Keldysh path integral associated to the Lindblad 
quantum master equation. The nature of the new FP 
precludes the use of more conventional classical dynam¬ 
ical field theories pioneered by Hohenberg and Halperin 
[35], and fully necessitates our quantum dynamical field 
theory approach. 

A platform for non-equilibrium quantum criticality - 
The starting point for our RG program is the quantum 
master equation governing the evolution of the density 
operator, p, of a one-dimensional (d = 1 ) bosonic field, 

$(x), 


d t p=-i[H,p\+£\p\. (1) 

Two-body collisions of strength A, among bosons (of 
mass m), are encoded in the Hamiltonian, H , while 
the Liouvillian can be decomposed into the sum of four 
dissipative channels C = ^2 a C a (a = P,l,t,d), with 
Ca[p\ = 7 a f x (L a (x)0 a (x) - ±{Ll(x)L a (x),p}), where 
local Lindblad operators incoherently create (destroy) 
single particles L p (x) = ffl(x) ( Li(x ) = 4>(x)), respec¬ 
tively with rates ( 7 ;), or destroy two particles, L t (x) = 
(f>(x) 2 , with rate 7 1 . The key element of our analysis is 
the Lindblad operator Ld(x) = d x cf>(x), which is respon¬ 
sible for single particle diffusion with rate 7 and which 
can be realised through microcavity arrays [9, 10] as por¬ 
trayed in Fig. 1 (see also [36]). 


In a simple mean field description, when the gain of 
single particles is balanced by two-body losses, a conden¬ 
sate cfo = (i4>{x )) with spontaneously chosen global phase 
can emerge [38-40], A crucial ingredient is now the ab¬ 
sence of particle number conservation due to pump and 
loss processes. For this reason, there is no sound mode 
with dispersion to ~ |q| (z = 1 ), and, instead, the canoni¬ 
cal dynamic exponent is z = 2. The effective phase space 
dimension is then D = d+z = 3, allowing for a condensa¬ 
tion transition in one dimension even when fluctuations 
beyond mean field are taken into account. 

A parameter regime of strong quantum diffusion can, 
indeed, disclose a quantum critical behaviour analogous 
to zero temperature quantum criticality, as we are go¬ 
ing to glean in the following, recasting the non-unitary 
quantum evolution encoded in Eq. (1) into an equiva¬ 
lent Keldysh functional integral formulation of dynamics 
[41, 42]. The quadratic part of the action, occurring in 
the Keldysh partition function, reads 

Skin = ^ (P c ,P q ) (pfl pTj (J) , (2) 

where 4> c and 4> q are the so-called classical and quantum 
fields, defined by the symmetric and antisymmetric com¬ 
binations of the fields on the forward and backward parts 
of the Keldysh contour [41]. In Eq. (2), P R = (P A y = 
id t + ( Kr — iKj)d 2 + ix is the retarded (advanced) in¬ 
verse Green’s function, while P K = i[f) — 2q ddf) is the 
Keldysh inverse Green’s function. In Eq. (2) we rela¬ 
belled the parameters in view of RG applications: at the 
microscopic scale, kjjv (the ultra-violet scale where our 
RG starts), they are expressed in terms of the couplings 
entering Eq. (1), K R \ kuv = Ki\ kuv = id\k uv = 7 d, 
X\k vv = and j\ kuv = 7 P + 7 ^ The existence of 

two independent Green’s functions, G R ! A and G K - an 
exclusive aspect of NEQ dynamics [41, 43] allows for 
a distinction between a ’’retarded mass”, x, which con¬ 
trols the distance from the condensation transition, and 
a ’’Keldysh mass”, 7 , which will play in the following the 
role of a temperature and which microscopically corre¬ 
sponds to a constant Markovian noise level induced by 
the environment. 

The role of canonical scaling of a non-Markovian 
quantum noise, P r (uj) ~ |w| (responsible for zero- 
temperature bosonic quantum phase transitions [41]), 
can be taken by the Markovian diffusive driving, 
Pneqid) ~ ^Adq 2 , hr a model with dynamical critical ex¬ 
ponent z = 2 (w ~ q z ). In our system, such quantum 
scaling regime and its associated NEQ fixed point, ap¬ 
pear then in the simultaneous limit 7 ; —> 7 P , and 7 1 , 
7 p —> 0, where diffusion becomes dominant over 7 in P K , 
and the mass gap closes (x -A 0). This double fine tuning 
is analogous to the simultaneous tuning of spectral gap 
and temperature to zero at equilibrium quantum critical 
points [25], and it opens the door to the realization of a 
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driven analogue of quantum criticality in our system. In 
passing we mention that the gapless nature of the NEQ 
drive PneqiQ —> 0 ) —► 0 , is due to the existence of a many- 
body dark state [ 11 , 12 ] - a mode decoupled from noise, 
located at q = 0 in our case. 


In a strongly diffusive, near critical regime, spectral 
and Keldysh components of the Gaussian action, Eq. (2), 
scale then alike, P R / A / K (q) ~ q 2 (in compact notation 
ypR/A/ A] _ 2 ). This fixes the canonical classical and 
quantum field dimensions to [<(> c ] = {(f> q \ = | and sets the 
canonical scaling dimension of quartic couplings to 2 — d 
(upper critical dimension, d c = 2 ). 


We now point out the two key scales which de¬ 
limit the scaling regime introduced here and the novel 
NEQ quantum critical point discussed later. The first 
can be gleaned from an analogy with equilibrium: In 
the quantum-classical crossover at finite temperature, 
the quantum scaling associated to equilibrium quantum 
phase transitions, persists at scales smaller than the de 
Broglie length, Lub. ~ [25], where temperature, 

T, cuts off coherent quantum fluctuations. Analogously, 
the novel NEQ quantum critical regime is delimited at 
low momenta by the Markov momentum scale, A m - the 
threshold where constant Markovian noise prevails over 
single particle diffusion, spoiling the diffusive scaling of 
the noise component of the quadratic action ( P K ). For 
A m we find the upper bound A m < 0.2Ag (see Supple¬ 
mental Material) and at distances larger than A~^, criti¬ 
cal properties are governed by a FP in the Kardar-Parisi- 
Zhang universality class [44] . The second key scale is the 
Ginzburg momentum scale, Aq — "ft/Ad- at momenta 
lower than this scale corrections to canonical scaling be¬ 
come effective, indicating the breakdown of a mean-field 
description [45]. According to this analysis, the novel 
critical behaviour manifests then in the momenta win¬ 
dow, A m < q < Aq. 

Non-Equilibrium Functional Renormalization- We 
now aim at determining the universality class, i.e. the 
full set of critical exponents associated to the quantum 
NEQ critical regime, which is technically characterized 
by the so-called Wilson-Fisher FP of the RG equations 
[45]. To this end, we dress the microscopic coefficients 
of Eq. (1) with RG corrections, employing a functional 
RG (FRG) suited for open NEQ quantum many body 
systems [21] (and previously developed for NEQ closed 
settings [46-48]). FRG allows us to interpolate from 
the microscopic dissipative action to the infrared effec¬ 
tive action, introducing an infrared regulator, R k , which 
suppresses stepwise fluctuations with momenta less than 
an infrared cutoff scale k. In this way we can approach 
smoothly the critical point where infrared divergences 
govern the physics. The FRG flow is based on a func¬ 
tional differential equation [32] for the effective action T*,, 


d k T k = | Tr 


(rf+i^)- 1 ^ 


, where the trace opera¬ 


tion, Tr, denotes summation over internal degrees of free¬ 


dom as well as summation over frequencies and momenta, 
and Tl the second functional derivative of the effective 
action with respect to the fields. In order to convert the 
functional differential equation for T*, into a closed set 
of non-linear differential equations for the RG running of 
the couplings (the beta functions [32, 45]), we provide a 
functional ansatz for Tfe = Sq^ = S k i n + Sint , where we 
systematically take into account in Si nt = Sh + S a all 
operators which are classified relevant according to the 
quantum power counting discussed above: 


(3) 

Sh and S a are respectively the hermitian and anti- 
hermitian parts of the interaction action. The potentials, 
U c = \u c {(j)* c (j) c — po) 2 and U q = ^u q {4> q 4 > q ) 2 , have associ¬ 
ated complex couplings u c ^ q = A Cj9 + iR c>q . which micro¬ 
scopically coincide with the parameters entering the mas¬ 
ter equation (A c | kuv = \\k uv = A, R c \k uv = R q \k uv = 
7 1 , 9i\k uv = 2 7 *)• The couplings g 2 and g 3 = X 3 + in 3 
are, instead, only generated in the course of renormal¬ 
ization. In Eq. (3) we introduced the condensate den¬ 
sity (resulting from balance of particles gain/losses), po, 
since in our practical calculations we approach the tran¬ 
sition from the ordered phase, taking the limit of the 
stationary state condensate po = (j)* c (j)c\ss = —i► 0 . 

In this way, we capture two-loop effects [32] necessary to 
compute the full set of critical exponents and thus deter¬ 
mine the universality class. We also rewrite the inverse 
R/A propagators allowing for a complex wave-function 
renormalization coefficient Z, P R = iZ*dt + K*d 2 , with 
K = Kr + iKj and Z = Zr + iZj, [49], whose anoma¬ 
lous dimension can acquire a real and imaginary part, 
r/z = Vzr + ipzr = —dtZ/Z. We mark that the quan¬ 
tum dynamical field theory, Sq^, has a richer RG op¬ 
erator content than a conventional equilibrium Martin- 
Siggia-Rose action [41, 43] or the semiclassical model in 
three dimensions for driven-dissipative condensation [ 21 ] 
(where 7 d = 0 , and u q = g\ = g 2 = g 3 = 0 from the 
outset - on the basis of their RG irrelevance). 

Rescaling all the couplings, {< 7 }, of Sq >k by the quan¬ 
tum canonical power counting, we find a FP solution, 
{ 5 *}, of the FRG beta functions in terms of the rescaled 
variables, {<?} (see [50]). The analysis in the vicinity of 
the FP gives access to the full set of critical exponents. 
We will use as a benchmark for the salient physical fea¬ 
tures of the quantum FP in d = 1 dimensions {D = 3), 
its semiclassical driven Markovian counterpart in d = 3 
dimensions [ 21 ]. 

A Non-Equilibrium Quantum Universality Class. 


S h = 
S a = 


/ x,t 


d u c - 

o 7 Pq 

d&c 


d K T* 

+ Wf’ 


d(j)q 


: + 


dCl 
d& 




- 4 [S 3 ( 0 <j) 2 - dliMq) 2 }- 
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Crit. Exps. 

V 

Vk r 

Vki 

Vzr 

Vzi 

h'ld. 

V-y 

Quantum 

0.405 

-0.025 

-0.025 

0.08 

0.04 

-0.26 

X 

Semi-classical 

0.72 

-0.22 

-0.12 

0.16 

0 

X 

-0.16 


TABLE I. Comparison between the critical exponents of the 
quantum and semi-classical driven dissipative models (taken 
from Ref. [21]). In the semi-classical scaling, 7 ~ k°, and the 
Markovian noise can acquire an anomalous dimension, 7 7 . 

(i) Non-Equilibrium Fixed Point- The key new prop¬ 
erty of the quantum FP is its mixed nature with coexis¬ 
tent coherent and dissipative processes, as shown in Fig. 
2. This new FP is less stable than the finite temperature 
FP, or the semi-classical one, since the additional fine 
tuning of the Markovian noise level is necessary to reach 
the quantum FP - as discussed above. 

In the domain of equilibrium phase transitions, the uni¬ 
versality classes of d-dimensional critical quantum sys¬ 
tems and of their classical d + z dimensional counter¬ 
parts [25, 51, 52] coincide. Table I compares the full set 
of critical exponents of the quantum transition and of 
its semiclassical driven-dissipative counterpart [ 21 ], elu¬ 
cidating that the analogy does not hold in the case of 
our NEQ setting. In the vicinity of the transition, the 
exponent (v) controlling the divergence of the correla¬ 
tion length of the Bose field, exhibits, for instance, the 
mismatch among the two critical characters. 

(ii) Absence of asymptotic decoherence- Persistence of 

quantum mechanical facets at criticality (for length scales 
shorter than is a common feature between our 

FP and equilibrium quantum critical points [25, 51, 52]. 
The low energy anomalous dispersion relation of critical 
modes, W& ~ k 2 ~ r,K i(c\ — m 2 ), encodes coherent effects 
(ci ,2 are two positive constants depending on the quan¬ 
tum FP), in contrast to the purely diffusive leading be¬ 
haviour of Wfc in the vicinity of the dissipative FP of the 
semiclassical model, u>k ~ —ik 2 ^ r,K i [21]. From an RG 
point of view, the exponent degeneracy, t]k r = r/Ry = 
—0.025, allows for a finite ratio of coherent propagation 
(Kr) versus diffusion (AT/), r = ~ k~ VK R +riK i , which 

is thus fully consistent with the results of Fig. 2, and in 
particular indicates absence of decoherence at long dis¬ 
tances. 

(Hi) Absence of asymptotic thermalization- The per¬ 
sistence of NEQ character at macroscales and the associ¬ 
ated non-thermal character of the distribution function, 
constitute the strongest evidence that the quantum uni¬ 
versality class found in this Letter cannot be related to 
its semiclassical driven Markovian counterpart in d + z 
dimensions, or to an equilibrium FP. 

To see this point, we note that the fluctuation- 
dissipation relation demands, in the 3-d driven- 
dissipative model, that the effective temperature Tq = 
\Z\*y - extracted from the infrared bosonic distribution 
function Fc(pJ,k) ~ —, is scale-invariant. This ex¬ 


Non Equilibrium Quantum Fixed Point 



FIG. 2. (Color online) In the quantum problem, all the 
rescaled couplings (we portayed some of the {g}) keep a non¬ 
vanishing real part at the FP and their RG flow freezes in the 
complex plane (as indicated by the red cross on the curved 
arrows). In the semi-classical problem, instead, decoherence 
forces asymptotically all the couplings to flow onto the imag¬ 
inary axis, and dynamics at infrared scales becomes purely 
dissipative. 

presses the principle of detailed balance of thermal equi¬ 
librium states (invariance of temperature under the sys¬ 
tem partition) in an RG language [21, 53]. Such circum¬ 
stance occurs at the semiclassical FP via the emergent 
exponent degeneracy ?y 7 = —pzR (cf. Tab. I) - the sys¬ 
tem thermalizes asymptotically. 

In the same spirit, if thermalization were to ensue 
close to the quantum FP, scale-invariance of the low- 
frequency distribution function, Fq(uj,k) ~ TQ J fc) (1 + 
7 */ 2 ) (Tq(A;) = \Z\"/dk 2 ), must be expected as a neces¬ 
sary condition. Specifically, replacing the bare scaling of 
the frequency ui ~ k z in Fq(oj, k), insensitivity to sys¬ 
tem’s partition would manifest in the exact scaling re¬ 
lation F q ~ k°. The absence of exponent degeneracy, 
r\ ld —t]zr (cf. Tab. I), signals scaling violation in the 
infrared behaviour of Fq ~ k~(vi d +rizR) ^ anc j accordingly 
the absence of infrared thermalization at the quantum 
FP. 

(iv) RG limit-cycle of Z- Finally, we notice that the 
peak of the spectral density - the imaginary part of 
the retarded single particle dynamical response, A(ui = 
Reo/fc) = R \z\P i m \j k ; is sensitive to oscillations present 
in Z ~ k~ VZR e ~ lr,zlt , induced by a non-vanishing rjzi 
(t = \og(k/Ac) is the RG flow parameter). Even if these 
RG limit-cycle oscillations occur with a huge period, , 
they are a remarkable signature of the novel critical be¬ 
haviour, since they prevent the possibility to have a real 
wave-function renormalization Z , contrary to what hap¬ 
pens for purely dissipative relaxational models [35] or for 
the semiclassical FP, where instead rjzi = 0 (cf. Tab I) 
[ 2 !]- 

Conclusions- We have shown that both quantum me- 
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chanical coherence and the microscopic driven nature 
of open quantum systems can persist close to a criti¬ 
cal point, in striking contrast with classical equilibrium 
and semiclassical NEQ fixed points. We have discussed 
the impact of this novel critical behaviour on the cor¬ 
relation length of the order parameter, on the distribu¬ 
tion function and on the spectral density. An important 
perspective direction is to study the effect of additional 
symmetries and conservation laws on the dynamical fine 
structure of this novel universal behavior. 
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RENORMALIZATION SCHEME 


The starting point for the renormalization scheme employed in this Letter is the Wetterich equation [32] 


dk^k 



Rk) 


(4) 


where the functional ansatz for Tk is the quantum NEQ Keldysh action, Sq^, of the main text. The trace operation, 

( 2 ) 

Tr, denotes summation over internal degrees of freedom as well as summation over frequencies and momenta; F), 
denotes the second functional derivative of the effective action with respect to classical and quantum fields, and 
finally Rk is the optimized cutoff [34], Rk = K(q 2 — k 2 )9(k 2 — q 2 ), where K denotes the bare kinetic coefficient. 


We now summarize the technical steps which lead us to disclose the Wilson-Fisher FP. 


We first introduce renormalized fields <p c = (j ) c , <p q = Z<f >„; this choice implies that the complex wave-function 
renormalization Z that multiplies the time derivative in p R 7 A , is absorbed in the field variables, and that kinetic 
coefficients and couplings are rescaled in the following way: K = K/Z , y <2 = 7 d/|Z| 2 , u c = u c /Z , u q = u q /{Z*\Z\ 2 ), 
9i = 9i/\Z\ 2 , 52 = <? 2 /l^| 4 > 93 = 93/Z 2 , 7 = t/\Z\ 2 . 

We project the functional flow encoded in Eq. (4), into the flow equations of the couplings, taking functional deriva¬ 
tives wrt to their linked quartic operators, and evaluating such derivatives on the classical and quantum background 
fields associated to condensation (ff> c = (f>* = (j} 0 and (j> q = <j>* = 0). 

The beta functions expressed in terms of the RG flow parameter, t = log (k/kuv)i read: 

<9 t A c = (d\ c = vzrK - Vzik c + AA c , 
d t K c = P Kc = qziK + VzrKc + Ak c , 
dt\ = 0\ q = 3 Vzr\ + VziKq + AX q , 

dfKq = /3 Kq = 3r]ZRKq ~ VZlX q + A Kq, 

dt.gi = p gi = ‘Z’nzRgi + A 51 , (5) 

dt92 = Pg 2 = dr] Z R92 + A 52 , 

d t X 3 = P\ 3 = ^VzrXs — 2ijzi r '3 + AA 3 , 
d t n 3 = /3 k 3 = 2r] Z RH3 + 2.q Z iX 3 + A k 3 , 
d t "1 = = 2r]zR r ) + Ay, 

where tjzr and rjzi are computed along the lines of Ref. [54] and they appear in Eqs. (5) as a consequence of the 
rescaling of the couplings through the complex wave-function renormalization coefficient, Z. The computation of 
AA C , ..., Ay has been carried extending the scheme of [54] to extract loop corrections of vertices which are quadratic, 
cubic, and quartic in the quantum fields. Further technical details of this procedure will be provided elsewhere [55]. 
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The NEQ critical properties under inspection are described by a quantum scaling solution of the flow equations, Eq. 
(5). We define three dimensionless ratios of coherent to dissipative couplings, r = Kr/Kj, rjj = X c /k c , r® = \ q /n q , 
and by virtue of the quantum canonical power counting, we introduce the dimensionless retarded tur = , 

and Keldysh mass, 7 = 7 ^ 2 , likewise the three dimensionless effective Keldysh masses fh R = 2 ^° , ttir = 2 ^2 , 
= 2 * d £° , and the dissipative rescaled couplings k c = ^ , k q = ^ ^, 32 = 37 • The flow equations of these 

dimensionless rescaled variables acquire a contribution from the running of the anomalous dimensions, r] a = —kdk In a, 
with a = Kr, Kj,"/d, which are functions of the rescaled variables too. 

Employing Eqs. (5), the flow of the ratios, rjj , r®, can then be written as 


d t ru = P ru = — (P\ c - ru/3 Kc ), 

K c 

d t r U = Pr? r = ttO 3 - r$P Kq ). 


( 6 ) 


For the ratio of real (Kr) and imaginary parts (Kj) of the kinetic coefficient, K , we have instead 

d t r = p r = r(rj Kl ~Vk r ), (7) 

where in the first equation we introduced the anomalous dimensions of the renormalized kinetic coefficient K, 


Vk r =Vk r 
riKj =Vk, ■ 


PZR 


Vzi 
r 

Vzr - Vzir, 


( 8 ) 


which are expressed in terms of tjzr , 3 zh and of the anomalous dimensions of the bare kinetic coefficients fjx R , 
fjK z - Vk r and fjx x are again extracted from the Wetterich equation, Eq. (4), following the procedure discussed in 
Appendix C of Ref. [54]. 


For the dimensionless masses of the problem, we have instead 


d t rh R = Pfn R = - (2 - t] Ki ) rriR + + -^p/3 Po , 

2 

dtmK = P^ = - (2 - r? 7d ) m l K + ^p(po/3 9l + 9iP Po ), 

2 

dtin\ = = - (2 - rj 7d ) m 2 K H- r^(poPx 3 + A 3 P Po ) 


(9) 


7 dk 2 

2 

d t m 3 K = p^ = - (2 - rj 7d )rh 3 K + -^~^{poPk 3 + k 3 /3 Po ), 
dtl = fa = - (2 - g ld ) 7 + ^ 2 /V 

where loop corrections to the condensate amplitude make their appearance in P Po (computed following Ref. [54]), as 
well the anomalous dimension for the diffusion coupling in the Keldysh sector 


d 9~id -fjZ R- 


( 10 ) 


r/~ /d is determined by loop corrections to 7 d, and it is computed generalizing the technique developed in [54] for the 
extraction of the bare kinetic coefficients fjx R , Vk t [55]. 


Finally, we need the flow equations of the following dimensionless dissipative couplings 

1 

dtRc = Pk c = ~ (1 _ 2 r]Kj + T 7H TS‘> I. Pk o > 


47t Kjk 


1 1 


dtii q — Pk q — (1 Vjd) K'q T ^ ’ 

dtg 2 =Pg 2 = -{ 1 + VKi - 2j/ 7 J 32 + -^^Tj,Pg 2 - 


( 11 ) 














Setting to zero the RHS of ( 6 ),(7), (9), (11), we find the quantum FP with coexistent coherent and dissipative pro¬ 
cesses: r* = 1.55, Vy = 0.397, r®* = -0.248, k* = 0.146, k* = -0.028 ,j% = 0.0047, rh* R = 0.900,7* = 0.023, m^* = 
0.22, m 2 ? = —0.0484, mg =0.803. 

Linearization around the fixed point provides the stability matrix, containing the information on the critical ex¬ 
ponents. The number of non-degenerate independent blocks in this matrix gives the number of independent critical 
exponents, which are five in our case (cf. Tab. 1 in the main text). 


ESTIMATE OF A M 


At equilibrium, A dB. is determined from the outset, since temperature is RG invariant. In our NEQ problem, the 
Markovian noise level j(k) is, instead, scale dependent and its value can be renormalized by loop corrections. 

After the RG flow enters the fluctuation dominated regime, k < Aq, loop corrections can sensitively affect the value 
of A m- The aim of the second part of our estimate is to extract the functional dependence of ”f(k) in proximity of the 
departure from the FP. In order to accomplish this goal, we consider the one-loop beta function for the Markovian 
noise level, where we rescale all the couplings besides 7 itself in view of replacing them by their FP values, 


kd k 7 


, 37 + 2fc 2 7 d 

91 3(1+ x *) 2 


( 12 ) 


where g 1 = 377 ^ and \ = K^k 2 • The solution of Eq. (12), 7 (k) = —b^dk 2 + Ck a , contains the FP-dependent 
parameter, a = 4 ^. 74)2 (with b = 3 ^+a) )i an< 3 th e integration constant C = (b + 7 *) 7 dAg +a , which lias 

been evaluated imposing the boundary condition for the differential equation Eq. (12) around the Ginzburg scale, 
7 {k ~ Aq) = 7*7dAg. 

Running couplings acquire rapidly their fixed point values when the flow enters the fluctuation-dominated regime, 
k < Aq. This justifies both the replacement by their FP values in Eq. (12), and the computation of the integration 
constant C. 

Solving 7 (k) ~ 27 dk 2 , we find then 


k — Aq 


/ 7 *+ 6 \^ 

V2 + b) 


(13) 


Replacing the FP values of the rescaled variables in Eq. (13), we find Am < 0.25Ag. 

We also performed a more detailed numerical estimate of Am using the above strategy but based on the two-loop 
beta functions. The result does not change significantly, Am < 0.2Ag, which is the value reported in the main text. 
For this reason, we presented here the technically less involved and physically more transparent estimate of A m from 
the one-loop beta function of 7 . 








